Sensitivity of ray dynamics in an underwater sound channel to vertical scale of 

longitudinal sound-speed variations 
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We investigate sound ray propagation in a range-dependent underwater acoustic waveguide. Our 
attention is focused on sensitivity of ray dynamics to the vertical structure of a sound-speed per- 
turbation induced by ocean internal waves. Two models of longitudinal sound-speed variations are 
considered: a periodic inhomogeneity and a stochastic one. It is found that vertical oscillations of a 
sound-speed perturbation can affect rays in a resonant manner. Such resonances give rise to chaos 
in certain regions of phase space. It is shown that stability of steep rays, being observed in exper- 
iments, is connected with suppression of resonances in the case of small-scale vertical sound-speed 
oscillations. 
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I. INTRODUCTION 

In recent two decades exponential divergence of ray 
trajectories with an infinitesimal difference in initial con- 
ditions, namely ray chaos iL>2a3jA,Mt stands as a top- 
ical problem in underwater acoustics. Practical sense of 
studying ray dynamics under chaotic conditions is mainly 
dealed with the hydroacoustical tomographv — monitor- 
ing of the ocean using sound signals [a, [3, [a, S ITU - 
Ray chaos raises difficulties in extracting an information 
about environment from data of long-range acoustical ex- 
periments using the traditional schemes of the tomogra- 
phy. Extreme sensitivity of chaotic rays to initial con- 
ditions sets the "predictability horizon" on ray motion. 
Furthermore, ray chaos leads to multiplying of eigen- 
rays and causes breakdown of the ray perturbation the- 
ory at long ranges fl^ IT^ . Nowadays development 
of efficient tomographical methods seeks the detailed de- 
scription of ray chaos and its manifestations in wavefield 
structure with the goal to obviate posed restrictions, if 
it is possible. The helpful feature in this way is sur- 
prising stability of the early portion of a received pulse. 
This phenomenon was observed in a number of experi- 
ments realized by NPAL group (see and references 
therein), for instance, in the well-known AET experiment 
[T^ . The numerical simulation of ray propagation in the 
AET environment shows that characteristic Lyapunov 
exponents of steep rays are less significantly than those 
of near-axial ones [ig. Such a picture indicates sup- 
pressing of diffusion in the range of large values of the 
action variable and strong chaos in the range of small 
ones. However, if the internal wave field is modeled as 
a sum of horizontal modes only, the range small values 
of the action corresponds to stable motion f3| . The pri- 
mary goal of the present work is to examine a role of 
vertical sound-speed variations in ray properties. The 



most attention is paid to the model of a waveguide with 
a deterministic periodic perturbation of a sound-speed 
profile. Although this model is idealistic, it yields a sim- 
ple description of chaotic dynamics in the framework of 
the KAM (Kolmogorov-Arnold-Moser) theory. The case 
of the stochastic perturbation is also considered to make 
our research more complete. 

The paper is organized as follows. In Section II we de- 
scribe briefly the model of a sound-speed profile, which 
we used in numerical experiments. Section III is devoted 
to analysis of ray dynamics with a periodic spatial inho- 
mogeneity of a waveguide. The sensitivity of ray prop- 
agation to the vertical scale of a perturbation is exam- 
ined by computing Poincare sections and timefronts of a 
received pulse. Section IV is concerned with ray propa- 
gation in a stochastically inhomogeneous waveguide. In 
Section V we summarize and discuss our results. The 
Appendix contains the exact expressions for the action- 
angle variables which was introduced for the model of an 
unperturbed waveguide. 



II. MODEL OF A WAVEGUIDE 

Consider a two-dimensional acoustic waveguide in the 
deep ocean with the sound speed c being a smooth func- 
tion of the depth z and the range r. One-way sound ray 
trajectories satisfy the canonical Hamilton equations 



dz _ dH 
dr dp ' 

with the Hamiltonian 



dp 
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dz ' 
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where n = co/c(z,r) is the refractive index, cq is some 
reference sound speed, p = nsini/) is so-called ray mo- 
mentum, and is ray grazing angle. As sound speed 
changes slightly with depth, i. e. \c(z,t) — cqI/cq <C 1, 
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one permits to rewrite the Hamiltonian (EJ in the parax- 
ial (small-angle) approximation 



Ac(z) (5c(z, r) 



Co 



Co 



(3) 



where Ac(z) = c(z) — cq, Sc{z, r) is a sound-speed pertur- 
bation along a waveguide, (Jcmax Acmax- It is clearly 
seen from lO that ray propagation in a range-dependent 
underwater sound channel is equivalent to motion of a 
nonlinear oscillator driven by an external force. Accord- 
ing to this analogy, a sound-speed profile plays role of a 
potential well, p is the analog to mechanical momentum, 
and r is the timelike variable. 

In Ref. 0, 0, ll^l we have introduced the realistic 
model of a background sound-speed profile 



c(z) = Co 



—az\f —az \ 

e )(e -7) 



0< z< h, 



(4) 

where 7 = exp(— aft,), /i = 1.078, a — 0.5 km ^, b — 
0.557, h = 4.0 km is the depth of the ocean bottom, 
Co = c{h) = 1535 m/s. 



waveguide over long distances and should be eliminated. 
Consequently, we can define the phase trajectory of the 
ray with E — as the "separatrix" of unperturbed mo- 
tion. As energy E isn't an invariant along the trajectory 
in a range-dependent waveguide, rays can transfer from 
one class to another. Whenever E becomes larger than 
zero, the ray reaches the bottom and escapes from the 
waveguide. Ray escaping is closely connected to trans- 
port properties of phase space [21I l23 | . In certain cases 
the portion of escaping rays enables to evaluate the rate 
of diffusion in phase space jlDl- In addition to that, we 
note that ray escaping takes place in waveguides having 
various physical nature, for instance, a similar effect is 
observed for the lasing modes in the chaotic dielectric 
microcavities 

For the unperturbed problem we are able to to reduce 
the ray equations to a more convenient form in terms of 
canonical action-angle variables The respective 

transformation is the following 



pdz, 1} = — I pdz. (6) 
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FIG. 1: Sound-speed profiles at the range R = 0: a) without 
perturbation, b) = 27r/1.0 km~^, c) fcz = 27r/0.2 km^^. 



For convenience we introduce the quantity E = 1 + 
H{e = 0) which we assign as the "energy" of ray oscilla- 
tions in a waveguide. Rays, propagating in the waveg- 
uide ^ over long distances, can be partitioned into two 
classes according to the value of E. The first class corre- 
sponds to those rays which don't reach neither the sur- 
face nor the bottom. For these rays E < Er, where Er 
is given by 



^r = y(l-A*)(l-7)- 



(5) 



Rays, propagating with reflections from the surface but 
without an interaction with the lossy bottom, relate to 
the second class. In this case energy satisfies to the in- 
equality Er < E < 0. Rays with E > reach the ocean 
bottom. Just as we consider the long-range sound propa- 
gation, number of the reflections from the bottom is large 
enough and, therefore, these rays cannot survive in the 



The exact analytical expressions for the action-angle 
variables and for the inverse transition to the variables 
(p-z) are presented in the Appendix. The action vari- 
able / measures the steepness of the ray and is equal to 
for the purely horizontal ray. Besides that, the action 
determines the mode the ray corresponds to by means 
the Bohr-Zommerfeld quantization rules 
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0<E<Er, 



Er < E <0, 



(7) 



where ko = 2nfl/co is wavenumber of the propagating 
wave, n is the career frequency. The angle variable i} 
serves as a phase of ray oscillations in a waveguide. 

The ray equations in terms of the action-angle vari- 
ables are written as follows 



dl _ 1 d{Sc) 
dr cq d'd ' 
dd _2ti 1 d{5c) 
dr D Cq dl ' 



(8) 



where the sound-speed perturbation Sc is expressed as a 
function of / and i9, and D is the ray cycle length (ray 
double loop). The ray cycle length in the waveguide © 
is given by the formulas 



27r 



E <Er 



D = I 



. A* + 7 - 2/i7 + 4^/6^ 
■K + Z arcsm 



E > Er. 

(9) 
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As it is shown from Fig. El D increases with increasing E 
from 37 km to 55 km Since the derivative dD/dE 
has two isolated zeroes near E = Er, the system of non- 
hnear coupled equations (QJ is locally degenerate. Such 
systems are known to exhibit chaotic behavior even un- 
der infinitesimal perturbations @. 




FIG. 2: The cycle length of the ray path D versus the energy 
E of ray oscillations in the unperturbed waveguide. 

In accordance with Fermat's principle, ray arrival time 
to a point R along a waveguide is calculated with the help 
of the Lagrangian L 

R R 

t=— I Ldr^— [{p^-H)dr. (10) 
Co J Co J 



At sufficiently long ranges, R/D 3> 1, the Lagrangian L 
may be considered as a function of the action 



L{I)=27T^-Ho{I). 



(11) 



Following to ray arrival time to the point R along 
a range-dependent waveguide is given by 



Co 

where (...) means averaging over r. 



(12) 



III. PERIODIC SOUND-SPEED 
PERTURBATION 

A. Theoretical analysis 

In this section we consider the periodic sound-speed 
perturbation given in the form 



dc{z, r) 



Co 



eV{z) cos {krV + kzz), 



(13) 



where e <IC 1 and V{z) is a smooth function playing role 
of an "envelope" of sound-speed variations. In numerical 



experiments we use the following expression for V{z) 



(14) 



where B = 1 km. The sound-speed perturbation II13|I in 
terms of the action-angle variables reads 



Sc{z, r) 
Co 



eV{I, z9)cos$ 



(15) 



where the phase $ — krr + kzZ varies along the ray path 
with the angle-dependent spatial frequency 



dr 



k{I, ^)=kr+ k,p{I, d). 



(16) 



The function V{I, d) can be decomposed into the Fourier 
series over the cyclic variable d 



V{I, ^)=Y1 ^"(^) cos mi?. 

m— 1 

Then the ray equations ^ take form 

m J 



rf7=-+2 EdZ^^^^^^ + ^^^^-'^- 



llT) 



(18) 



(19) 



dz \ 

- fc. 51 ^"^ (^in *m + sin *m) j , 

where u) — I-k jD and — m-d ± krr ± kzz{I, t9). The 
stationary phase conditions d'i'^Jdr ~ imply the first- 
order nonlinear resonances 



muj{I) - kr - kzp{I, d) ~ 0, 



muj{I) + kr + kzp{I, 1?) ~ 0. 



(20) 



(21) 



We intent to investigate how interrelation between dif- 
ferent length scales of a perturbation influences the rays. 
Three qualitatively different regimes of ray motion can 
be separated according to the ratio of kz and kr : 

1) l^zPrnaxI kr', 

2) I^^J'maxI — kr, 

3) IfczPmaxI 3> kr, 

where Pmax is the ray momentum at the channel axis. 

The first regime corresponds to a horizontal sound- 
speed perturbation (see, for instance 0,SE,l23|)- In 
this case the resonant condition is the simplest 



muj{l — /res) = Ikr, 



(22) 
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where / and m are integers. The action of a ray under 
a resonance ll22|l oscillates with range. These oscillations 
are described by the universal Hamiltonian of nonhnear 
ill 



resonance 



Hu^m{- \uj'j{I,^,)\ {Aiy + e|V;„| cos^'; 



(23) 



where u;j(/res) = d^{Ives) / dl , A/ — I — I^es- The width 
of an isolated resonance in terms of the spatial frequency 
of a ray trajectory is approximately estimated as 



-■/I A/jiiax 



(24) 



where A/„iax is the width of the resonance in terms of the 
action variable. In accordance with lfT2|l . if the distance 
between the source and the receiver is large compared 
with the period of action oscillations, all the rays, being 
captured in a given resonance, have close arrival times. 
Consequently, sharp peaks and gaps arise in distribution 
of ray arrival times at a receiver 

In accordance with Chirikov's criterion (2^, global 
chaos arises if 



Suj 



1. 



(25) 



i. e., if two nonlinear resonances, centered at uj and lu + 
Suj, overlap. Those resonances that overlap sHghtly form 
islands in phase space, areas of stable ray motion in the 
chaotic sea. The distance between the resonances of the 
TO-th and m + 1-th orders in terms of spatial frequency 
is equal to 



k k 

OUJ = 

771 m + 



UJ 

— ocD- 
k 



(26) 



and decreases rapidly with increasing D. Hence chaotic 
properties of rays under a horizontal sound-speed pertur- 
bation depend on the shape of the unperturbed profile 
defining the function D{E) [28j. 

Under the second regime the condition Ij21ll can be 
satisfied, that corresponds to capturing a ray into the 
resonance with vertical oscillations of a sound-speed per- 
turbation. This phenomenon is similar to the "wave- 
particle" resonance [2fll l8n| . Hereafter we will call this 
sort of resonances as "vertical resonance". Vertical reso- 
nance arises at some fixed point where the ray momen- 
tum is equal to 



muj{I) 



(27) 



Define the /^-vicinity of a given ray as a small region 
along a ray trajectory, where the detuning of vertical 
resonance doesn't exceed some small value fi, i. e. 



lAI = 



(28) 



Influence of vertical resonance depends on how long the 
ray passes through the ^-vicinity. Consequently, the 



maximal influence is expected, when the /x-vicinity in- 
cludes the channel axis where dp/dr = and p = max. 
Therefore, vertical resonance affects mainly those rays 
which obey to the equation 



Pn 



muj{I) 



(29) 



If UJ is small compared with kr and one takes into ac- 
count that the Fourier amplitudes Vm decay rapidly with 
growth of 777, the resonant condition (I29II has the simplest 
form 



Pn 



(30) 



Interplay of vertical and usual "horizontal" (i. e. relating 
to the condition lO) resonances distorts resonant tores 
[2^. Such a distortion can enforce close resonances to 
overlap, that causes emergence of global chaos. 

Under the third regime behavior of near-axial and 
steep rays is qualitatively different. A perturbation 
causes strong distortion of the sound-speed profile near 
the channel axis, and the derivative duj/dl vanishes in 
the range of small values of the action, that is a sufficient 
condition of chaos of near-axial rays Furthermore, 
fast depth oscillations of a perturbation lead to occur- 
rence of microchannels near the waveguide axis, and in- 
stability is amplified by irregular jumps of rays from one 
microchannel to another. 

On the other hand, fast oscillations of a perturbation 
along steep rays can be eliminated using the averaging 
technique and, therefore, the problem can be reduced to 
the integrable one 

^liaHi- The averaging technique 
the spatial frequency u to be small compared with 
Certainly, this condition doesn't hold near the refractive 
ray turning points where fc = (under reflections from 
the surface \p\ ^ and k can be large enough). In these 
regions influence of horizontal resonance II22II has to be 
taken into account. Nevertheless, the size of the "reso- 
nant" region, being estimated as Ap ~ k~^ , is small for 
large k^. Moreover, it can be shown that the impact of 
horizontal resonance ceases with increasing kz. Let us 
rewrite the range-dependent term in the Hamiltonian in 
the form 

Hi = V{1, 1?) exp[ikzz{1, 7?) + ikrr] + c. c. (31) 
Then we get Fourier-expansion of the function 2(7?) 



z(7?) = Zq + ^ Zn cos 777?, 



(32) 



and substitute it into the equation (j^lll . Leaving the 
superior term with 77 = 1 only, we can represent li;^l|l in 
the form 

Hi — Veff exp {ikrr) + c. c, (33) 

where Veff is expressed in terms of the Bessel expansion 

00 

Ve{{^V{I,d)eM^kzZo) i'Ji{kzZi)eMii^)- (34) 

I — — 00 
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According to H24II and the resonant response of a 

ray can be expressed as 



In the limit kzZi 
asymptotics 



AI ^ VKff. (35) 
oo the Bessel functions have the 



JiikzZi) 



1 



(36) 



Therefore, the resonant response AI decreases with in- 

; ; -1/4 
creasing kz as Kz 

Thus the criterion of appHcability of the averaging 
technique can be formulated in the following way: in- 
equalities 



dV 



dz 



<t:\kzV\, uj<t:kz\p\, 



(37) 



must be satisfied for that range of the depth, where the 
function V{z) differ considerably from 0. The first in- 
equality in ll37ll infers insignificance of resonances along 
the trajectory. The second one requires the ray momen- 
tum to be large enough. The averaging technique reduces 
the ray equation written as 



dr'^ 



1 

Co 



d{Ac) d{Sc) 
dz 



into the "smoothed" form 



dr"^ 



1 

Co 



d(Ac) 
dz 



M 



dz 



d{Sc) 
dz 



(38) 



(39) 



where M is some averaging operator. In the vicinity of 
the channel axis the spatial frequency k varies slowly and 
the averaged term in ll39|l can be expressed as follows 



M 



djSc) 
dz 



e^klV dV 
2fc2 



(40) 



B. Numerical simulation 

The Poincare map is known as an efficient method 
to separate stable and unstable rays in the underlying 
phase space. We constructed Poincare maps in the polar 
action-angle coordinates with different ratios of kz and 
kr ■ The result is demonstrated on Fig. kz varies from 
to 27r/0.2 km~^, while kr is chosen to be of 27r/10 km~^ 
for all computations. Other parameters of the perturba- 
tion have the following values: s = 0.0025 and B = 1 km. 

When kz is very small or equal to zero, ray dynamics 
is predominantly regular (Fig.|SK)- Chaos occurs only in 
the isolated chaotic layer near / = I{Er), where nonde- 
generacy violates. Nevertheless, inclusion of even slow 
vertical oscillations broaden the chaotic layer greatly, 
that is shown in Fig. Eh- The so-called "chaotic sea" 



emerges and chaotic diffusion in phase space becomes 
unbounded and global. The chaotic sea is rarely filled by 
points due to strong ray escaping. Such amplification of 
ray escaping is induced by vertical resonance II^OII which 
affects especially the steep rays satisfying the condition 



bmax(/)|-T^=0.2. 



(41) 



With further increasing kz the most steep rays exhibit the 
transition chaos - stability. It is shown in Fig. Eh, where 
the Poincare map for the case oi kz = 27r/1.0 km~^ is de- 
picted. Now the chaotic sea is separated into two parts, 
the internal layer and the external one (Fig. Et)- Both 
these layers are restricted by invariant curves. Complete 
stabilization of the steep rays takes place in the case of 
kz = 27r/0.2 km~i (Fig. 131). With that, the internal 
chaotic layer occupies entirely the range of small values 
of the action. 

In Fig. 01 fragments of three ray trajectories with kz, 
varying from 27r/2.0 km~^ to 27r/0.2 km~^, are shown. 
Initial conditions and other parameters of the sound 
channel are the same for all trajectories. All the rays 
are launched from the channel axis with the same start- 
ing momentum p{r = 0) = 0.03. The ray with kz — 
27r/2.0 km~^ is completely stable, that is in agreement 
with Fig. 03. In the case of kz = 27r/l km~^ the ray 
belongs to the internal chaotic layer and, therefore, is 
chaotic (this case corresponds to Fig.^fc), but its irregu- 
lar behavior is not so apparent, as that of the ray with 
kz = 27r/0.2 km~^. In the latter case the "elevating" 
effect takes place - the ray, being trapped into a mi- 
crochannel near the lower turning point, leaves it only 
near the upper turning point. The action of an "elevated" 
ray insreases appreciably. 

The qualitative description, being provided by the 
Poincare maps, enables to explore the properties of the 
wave field in the time - depth plain at a receiver, the so- 
called timefront. Timefronts were the main measurable 
characteristics in the recent long-range experiments fl^ . 

Chaotic dynamics of rays leads to divergence of their 
trajectories and, therefore, arrival times of chaotic rays 
with close launching angles can be very different. Besides 
that, topology of the Poincare maps (Fig. shows that 
only the rays with certain values of the action perform 
chaotic motion, whereas the rest part of phase space cor- 
responds to stable rays. Therefore, arrival time spreading 
is significant only in some finite segments of a timefront 
of a received pulse. 

We computed timefronts with different values of kz. 
Each timefront is modeled as superposition of arrivals 
of 50000 rays started from the channel axis. The ini- 
tial momentum p is uniformly distributed in the interval 
[-Pmax : Pmax], where pmax p{E = 0) ~ 0.22. The 
timefronts for kz — and kz — 27r/0.2km^^ are pre- 
sented in Fig. El In the case kz = the time spreading is 
weak. There are two prominent gaps and densely filled 
region amid them in arrival pattern, those correspond to 
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FIG. 4: Ray trajectories with different values of vertical 
wavcnumbcr of the perturbation k^. All the rays arc launched 
from the axis of the waveguide with one and the same starting 
momentum p{r = 0) = 0.03 



the horizontal resonance with I = 1 and rn = 5. Such 

gaps arc absent in the case of = 27r/0.2km^^, that 
is connected with suppressing of horizontal resonance in 
the regime of fast vertical oscillations of a sound-speed 
perturbation. In this case the final portion of the pulse is 
chaotic and unresolved, that is a consequence of strong 



chaos of near-axial rays with smah values of the action. 



IV. RANDOM SOUND-SPEED 
PERTURBATION 

Internal waves in the deep ocean are known to have 
continuous spectrum and must be considered rather as 
a stochastic process than a deterministic one. So the 
fohowing question arises: in what extent results of the 
previous section can be referred to acoustic propagation 
in a randomly perturbed sound channel? Our goal now is 
to investigate how ray dynamics depends on the vertical 
scale of random longitudinal sound-speed variations. 

In the present section we consider the sound-speed per- 
turbation given in the form 

(5c(^,r) = £Co-|e-2-/^C(^,r), (42) 

where a random function ^(z, r) is expressed as 

^{z, r) = cos [tt (ue-"' + \/2^4|(r))] |(r), (43) 

where A = 1, v = 1 km^^ and ^(r) is a random function 
with normahzed first and second moments 



0, 



(44) 
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FIG. 5: Timefront of a received pulse at the range 1000 km 
in the periodically perturbed waveguide: a) = 0, b) = 
27r/0.2km-^ 







a) 

■ 






FIG. 6: Sound-speed perturbation vs depth: a) it = 5, b) 
u = 20. 



A function ^(r) is modeled as a sum of 10000 randomly- 
phased harmonics with wavenumber distributed in the 
range from 27r/100 km~^ to 27r/l km~^. Spectral density 
of ^(r) decays with k as k~^. Single realizations of the 
perturbation with u = 5 and u — 20 are presented in 
Fig.H 

We integrated numerically the ray equations with the 
parameter u varying from 5 to 20. Figure represents 
the ray action at the range 1000 km as the function of the 
starting action Iq . In the case of slowly varying vertical 
structure with u = 5 (upper plot) chaoticity grows grad- 
ually with increasing Iq. Flat rays don't perform extreme 
sensitivity to initial conditions and dependence is contin- 
uous in the range of small Iq ■ Isolated continuous pieces 
belong to bundles formed by stable or weakly chaotic 
rays with close launching angles and converging travel 
times, namely coherent clusters Coherent clusteri- 

zation is a sort of cooperative effects in randomly-driven 
nonlinear Hamiltonian systems Creation of coher- 

ent clusters results from the resonant interaction of un- 
perturbed ray motion with low-frequency components of 
a random perturbation. In point of fact coherent clusters 
can be considered as counterparts of islands of stability. 
Continuous pieces melt with increasing /q and the de- 
pendence I{Io) become almost irregular, that attests ray 
divergence and chaos. 



In spite of the case u = 5 there are no any remarkable 
coherent clusters in the case of u = 20. The dependence 
/(/o) is chaotic for all values of the starting action. It 
should be noted that the region of small values of the 
action is more fuzzy than the region of large ones. More- 
over, spreading of values of the action depends slightly 
on Iq for steep rays, that implies hornogeneous chaotic 
diffusion. According to Virovlyansky Ho, 113' if ^^J dif- 
fusion in phase space is homogeneous and can be treated 
as a Wiener process, early arrivals are expected to be 
well-resolved. This statement is confirmed by the time- 
fronts computed with u = 5 and u = 20 (Fig. EJ. If 
u = 5 (upper plot) all the arrivals are well resolved, but 
their distribution is nonuniform: coherent clusters, look- 
ing as concentrations of plotting points with extremely 
small time spreading, alternate with rarely fihed regions. 
On the other hand, all the timefront branches are fuzzy 
in the case of u = 20. Time spreading, however, abates 
with decreasing the arrival time and is less than distance 
between neighboring branches in the early portion of the 
pulse. Locations of branches and the envelope of a time- 
front in the early portion are coincide with ones for the 
unperturbed problem. As long as these characteristics 
can be considered as "marks" of environment, those are 
relevant for the long-range acoustic thermometry [l^lSl 
or monitoring of the large-scale ocean variability • 
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FIG. 8: Timefront of a received pulse at the range 1000 km 
FIG. 7: Ray action at the range 1000 km as a function of its randomly perturbed waveguide: a) m = 5, b) m = 20. 

initial value in the randomly perturbed waveguide, a) ?i = 5, 
b) M = 20. 



SUMMARY AND DISCUSSION 



The stability of early portion of a received pulse can be 
clarified in the same way as that it has done in the case of 
the periodic perturbation in the Section III. Correlation 
length of sound-s^eed variations along a ray path is given 
by the formula (4(i | 



L 



1 -p2 



'1/2 



(45) 



where Lr and Lz are the horizontal and the vertical cor- 
relation length scales of sound-speed variations, respec- 
tively. If Lz is small enough, correlation length is ex- 
pressed as L = Lz/p and, therefore, is minimal for steep 
rays. Thus, if the criterion lIHTIl . being expressed as 



dV_ 

dz 



\p\ 

^\kzV\, CJ<^, 



(46) 



is satisfied, fast oscillations of the perturbation suppress 
each other and diffusion in phase space weakens, that re- 
veals itself in diminishing time spreads in the early por- 
tion of a received pulse. 



In the present paper we have investigated in what man- 
ner stability of rays in an underwater sound channel de- 
pends on the vertical scale of a sound-speed inhomogene- 
ity induced by internal waves. The main result is that 
fast depth oscillations of a sound-speed perturbation can 
lead to stability of the steep rays forming early portion 
of a received pulse in experiments. We found that a ray 
can be captured in resonance with vertical sound-speed 
variations, which we call as vertical resonance. If the 
length scale of vertical variations is large enough, verti- 
cal resonance enforces strong chaos and escaping of steep 
rays. 

Our results rely on the ray approximation without re- 
ferring to the pulse characteristics, but we suggest that 
they can be worthwhile for understanding some features 
of the propagation of a monochromatic acoustic wave. 
This suggestion is mainly based on the results presented 
in the recent paper there the task of the ray /wave 
method correspondence was studied. As was reported 
in improving of ray predictions requires to smooth 
those features of internal wave fine structure, which have 
vertical length scale less than the threshold of wavefield 
responsibility defined by the formula 



An 



Co 



f2 tan ( 



(47) 
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where is the career frequency and tan 0max is the maxi- 
mal grazing angle. For = 75 Hz this scale of roughly of 
110 m. On the other hand, we have found ray dynamics 
to be very sensitive to the fine structure of the internal 
wave field, that implies the arrival pattern to change with 
varying frequency. For instance, the early portion of a 
received pulse can lose stability at very low frequencies, 
when the threshold given by ()47ll is large. This seems 
to associate the anomalous low-frequency sound attenu- 
ation with strong ray escaping. 

However, there are many questions concerning struc- 
ture of late arrivals. First, Amin is large enough for flat 
rays, so they are expected to be insensitive to the small 
vertical features of the background sound-speed profile. 
Second, it is impossible to use geometrical acoustics to 
describe the propagation along the axis of a waveguide 
due to multiple caustics i4S|. On the other hand, forming 
of coherent ray clusters should explain the well-resolved 
coherent late arrivals which were observed at 28 Hz in the 
AST experiment 0|- the following question arises: 
in what extent complicated structure of the arrival final 
is linked with ray chaos? This topic will be aim of our 
future research. 



ACKNOWLEDGMENTS 

This work was supported by the project of Far East- 
ern Branch of Russian Academy of Sciences "Acousti- 
cal tomography at long ranges under conditions of ray 
chaos". We wish to thank V.A. Bulanov, S.V. Prants 
and A.O. Maksimov for helpful discussion during course 
of this research. 



the following formulas 



J ^ P{z ^ 0) ^h_l /^ + 7 _ M + 7 
TTfl a \ 4 27r 




(51) 



d = 



TT — iJr 



. M + 7 - (2/i7 - AE/b'^) 
arcsm 



Q 



7r-i9,. 



TT 

2 + 



. u + 7 - (2u7 - AE/b'^) 
arcsm — — - — 



, p<0. (52) 



In l|5Tll - lf52|l we used the notation 

> + 7- 2fj.'y + AE/b"^ 



(53) 



Under reflections, the ray momentum is given by the 
formula 



p{z - 0) - ^2E-b^{^l-l)iJ-l). 



(54) 



The inverse transformation for the rays, propagating 
without reflections from the surface, is expressed as fol- 
lows 



Appendix 



z(I, i}) = - \n 
a 



1 a^b^ (/i + 7 - Q cos I?) 
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(55) 



The action variable for the rays with E < Er is given 
by the formula 



2E 



(48) 



The respective expression for the angle variable is the 
following 

i9 = ±- =F arcsm — , (49) 

2 Q 

where the quantity Q is given by the formula 



8E_ 



(50) 



The upper and the lower signs in l(i9jl correspond to the 
case of p > and to the case p < 0, respectively. The 
action and the angle for surface-bounce rays are given by 



Lo Q sin 



a (/X + 7 — Qcosi?) ' 



where 



a6(M + 7)_„2,, 



a;'(/) 



(56) 



(57) 



The quantity uj' coincides with spatial frequency of ray 
oscillations uj for the rays propagating without reflections 
from the surface. Position and momentum of the surface- 
bounce rays are expressed as follows 

1 a^&nA^ + 7 + Qcos(^(^ + 7r))] 

- TT < I? < 0, 

1 a^b^[y^ + ^-Qcos{^d + dr)] 

z{I, d)^- In — 2 , 

a 2u;'^ 

< I? < TT. 

(58) 
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p{I, = 



a;'Qsin(^(^ + 7r)) 



/i + 7 + Qcos(^(79 + 7r)) 



TT < ?9 < 0, 



p{I, t?) 



(59) 



< 7? < TT. 
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